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SU(N C ) colour-singlet restriction, along with flavour and spin symmetry, on thermal partonic 
ensemble is shown to recombine the partons with internal colour structure into colour-singlet multi- 
quark clusters which can be identified with various hadronic modes at a given temperature. This 
can provide a possible basis for recombination model for hadronization of quark-gluon plasma. This 
also leads to a natural explanation for the ratio of (anti)protons to pions and the quark number 
scaling of the elliptic flow coefficient in relativistic heavy-ion collisions. 
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Introduction: By now it is generally believed that the 
deconfincment of strongly interacting matter has been 
achieved in relativistic heavy- ion collisions Remark- 
able confirmations for this premise come from elliptic flow 
of hadrons and its scaling with the number of valence 
quarks pL jet quenching and radiation of thermal 
photons 4j. However, the results from RHIC experi- 
ments also revealed some interesting facts that the nu- 
clear suppression factor depends on the hadron species Q 
and the proton to pion ratio [f| has a plateau around 
unity in the transverse momentum range (2 - 4) GeV/c. 
Recently, it has been argued in Refs. [y, 0] that the 
hadron production at low momenta in a dense medium 
takes place through the recombination of partons which 
explains some of the surprising results of RHIC experi- 
ments. 

Generally, hadronization by recombination is a "two 
to one" process in a medium where the (anti)quarks are 
effective degrees of freedom and gluons are dynamical 
ones that disappear at hadronization. This simply boils 
down to the fact of a correct counting of quantum states 
and momenta within a dynamical theory. However, it 
is an extremely difficult task within the dynamical QCD 
and various models [ESULSIlil have been formulated to 
describe the hadron productions in heavy-ion collisions. 

Recently, under a sudden approximation the recombi- 
nation is considered with a perturbative quark, i.e., mini- 
jet and a thermal (anti)quark Q whereas in Ref. [|| it is 
with thermal (anti)quarks. In Ref. Q it is argued that 
an additional contribution is necessary for the transverse 
momentum spectra of hadron at the transition region 
between the thermal recombination and the individual 
fragmentation. Such contribution comes from the recom- 
bination of a minijet and a thermal quark. In Ref. Q, 
on the other hand, it is argued that for momenta below 



5 GcV/c the thermal recombination dominates whereas 
beyond 5 GeV/c the fragmentation of independent mini- 
jet dominates the hadron production. Also, the com- 
petition between the recombination and fragmentation 
pushes the onset of fragmentation to relatively higher 
transverse momentum of 5 — 6 GcV/c. This indicates 
that the quark recombination phenomena for hadronisa- 
tion of quark-gluon plasma remain an open as well an 
interesting problem. 

In the present article, we show that the recombina- 
tion phenomena arise spontaneously upon application of 
colour-singlet projection operator on the partition func- 
tion for an assembly of quarks and antiquarks having 
internal colour structures. Such recombination phenom- 
ena naturally explain the baryon-to-meson ratio and the 
azimuthal anisotropy of hadron distributions that scales 
with the number of valence quarks. 

Quantum statistical mechanics and colour- singlet ensem- 
ble: The statistical behaviour of a quantum gas in ther- 
mal equilibrium is usually studied through an appropri- 
ate ensemble. In general one defines a density matrix for 
the system as 



P {f3) = exp(-/?£0 



(1) 



where f3 = l/T is the inverse of temperature and H is 
the Hamiltonian of the physical system. The correspond- 
ing partition function for a quantum gas having a finite 
volume can be written as 
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where |n) is a many-particle state in the Hilbert space 
Ti. and V is the projection operator for any desired con- 
figuration. We propose to consider the statistical ther- 
modynamical description of a quantum gas consisting of 
quarks and antiquarks, such that the underlying sym- 
metry amounts to reordering of the partition function in 
terms of the colour- singlet multi-quark modes at a given 
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temperature. We also assume here that the gluons are in 
thermal background and the partonic matter is mainly 
composed of quarks and antiquarks. 

Now for a symmetry group Q (compact Lie group) hav- 
ing unitary representation U (g) in a Hilbe rt sp ace H, the 
projection operator can be written as [HI, Elf 

V^dj f dn(g) x *(g)U(g) , (3) 
Jg 

where dj and \j are : respectively, the dimension and the 
character of the irreducible representation j of Q. dfi(g) 
is the normalised Haar measure in the group Q. The 
symmetry group associated with the colour-singlet con- 
figuration of the system is SU(N C ), N c is the number 
of colour corresponding to fundamental representation. 
For the SU(N C ) colour-singlet configuration dj = 1 and 
Xj = 1 ■ Now the colour-singlet partition function for the 
system becomes, 

Zs^Trlf dfi(g) U(g) cxp (-/?#) J , (4) 

\JSU(N C ) ) 

where U(g) can be thought of a local link variable that 
links the (anti)quarks in a given state of the physical sys- 
tem. Now, the trace in a Fock space results in a product 
of the fcrmionic determinant over momentum modes and 
colour degrees of freedom, which we shall see later. Thus, 
the interchange of the 'Tr' and the integration will lead 
to an overall product over momentum and color states 
and one can write the (0| as 

2s = TT / d ^) Tr (&(g)exp(-p6)) ■ (5) 

JSU(N C ) V 7 

We neglect the mutual interactions among the con- 
stituents, although they must interact in order to come 
to a thermal equilibrium. One can imagine a situation 
by first allowing them to come to a thermal equilibrium 
and then slowly turning off the interactions [Tj|. Such 
a simple thermodynamic description is often useful for 
various physical systems (e.g., electrons in metal, black- 
body photons in a heated cavity, phonons at low tem- 
perature, neutron matter in neutron stars, etc.). The 
full Hamiltonian is then the sum of the Hamiltonians for 
each species, i.e., quarks and antiquarks as H = H q + Hq 
with Hi = hi — HiNi, in the grand canonical ensemble 
with the usual meaning of \ii and iVj. 

Now, the Hilbert space TL of the composite system has 
a structure of a tensor product of the individual Fock 
spaces of quarks and antiquarks as TL — TL q ® Hq, where 
the subscripts q and q denote the quark and antiquark, 
respectively. Due to this the partition function in (J5]) 
in Hilbert space TL decomposes into the product of two 
traces [HI, H, EM Eil in their respective Fock spaces as 

Zs=\[ f d M (<7) Tr(u q (g)e-e*«)V:(u q -(g)e-^)(&) 

JSU(N C ) V 7 V 7 



For simplicity, we approximate the local link variable 
U(g) to be diagonal matrix related to only the diago- 
nal generators in colour space associated with maximal 
abelian sub-group (Cartan space) [ll| of SU(N C ). As for 
example, SU (3) colour gauge group has only two param- 
eter abelian sub-group associated with the two diagonal 
generators that would characterise U(g), including its di- 
agonalisation as can be seen below. Under this approxi- 
mation the Haar measure correspon ding to SU (N c ) can 
now be written in the Weyl reduced [ill [l2j form as 




where J is the Jacobian of transformation (also known as 
Vandermonde determinant [TtJ )■ Qi is a class parameter 
characterizing the group element g such that U(g) can 
be diagonalised. It also obeys the periodicity condition 
Sfci ^ = (mod27r), which ensures that the group ele- 
ment is SU(N C ). Thus, it is obvious that the product of 
two Fock space traces in © has to be a class function, 
which can be obtained below using the diagonalization 
condition in maximal abelian sub-space of SU(N C ). 

Now, in each Fock space there exists a basis that diago- 
nalizes both operators as long as Hi and Ui(g) commute. 
Let \a,a) be the one-particle states in such basis, where 
a labels the eigenvalues of Hi (including a possible de- 
generacy besides the one associated with the symmetry 
group Q) and a labels those of Ui(g). One can write the 
diagonalised cigenstatc as 

(a',a'\Ui(g)Ai\a,o-) = 

where Ai a and Ri aC r(g) are, respectively, the eigenvalues 
of Ai and U Aq) . Then following standard procedures one 
can obtain [lj, EH 

Tr (Uie-^ =exp [tr Q tr c In (l + R l {g)e~ l3A '-)] . (9) 

Note that 'tr a ' is the trace over the momentum state 
a whereas 'tr c ' is the trace over the colour degrees of 
freedom in the same momentum state a. In fundamen- 
tal representation tr c R(g k ) = exp(ik8i) along with 
R k (g) = R(q k ). This can be related to the Polyakov 
Loop 0, EH in Polyakov gauge as C = tr c (L)/N c = 
tr c R(g k ), where L is the thermal Wilson lines defined 
by the temporal gluons in Euclidian time. This corre- 
spondence is due to the choice of diagonal U(g), which 
resembles the Polyakov Loop matrix in Polyakov gauge, 
supplemented with the diagonalization condition in ([S]). 
Now one can write the product of two traces in ([6]) as 

^(u q {g)e-^)Tt{Uq{g)e-^)= exp(9), (10) 
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where it is easy to show that 

N c N s N f 

e = x)£££[ in ( i+eitf,e " /,(e,a "'*' ) 

a l s f 



(11) 



Here the flavour (/) and spin (s) summations are intro- 
duced, where Nf and N s are, respectively, the number of 
flavour and spin degrees of freedom. The equation (fTTj) 
clearly indicates the superposition of two Fock spaces 
[fcrmionic and antifcrmionic determinants] with parti- 
cles having same momentum and internal colour struc- 
ture obeying the quantum statistics. The momentum 
states which do not satisfy this are automatically elim- 
inated by the diagonalization condition in fSJ). This es- 
sentially amounts to a stacking of same momentum par- 
ticles. It is also evident that even if we had started 
by considering a free gas of quarks and antiquarks, the 
colour-singlet restriction related to the Polyakov Loop 
in Polyakov gauge links the effective degrees of freedom 
(anti)quarks with the surrounding thermal bath through 
the temporal gluons. As we will see, this allows an in- 
teresting scenario of recombinations. Now, the single- 
particle energy eigenvalues in a given state are same for 
q and q as Ci a = \J pf a + mf . However, their occupation 
energies in a given state differ by their chemical poten- 
tials as fig = —fj,q = —(J,. For convenience [l5| , we make 
a substitution £ = —ij3jJL and e qa = e a in (fTTj) which then 
becomes 



e = inJIII v 



a I 



V= le - ^" (2cosh/?e Q + 2cos(0j+f)) ,(12) 

where the determinant V is a class function. We would 
like to note that the class parameters in colour space get 
associated with the imaginary chemical potential £ due 
to the choice of diagonal U (g). As a consequence of which 
we will see below that the colour factor N c appears with 
the chemical potentials indicating the number of valence 
(anti)quarks. 

The partition function in can now be written as 



n 



(13) 



SU(N C 



The integrations on class parameters in Q13p arc now per- 
formed exactly by using the properties of the Jacobian 
and an orthonormal polynomial method [13] ■ The loga- 
rithm of the partition function [16( for N c = 3 and two 
massless quarks (Nf = 2) in the infinite volume, V, limit 
reads as 



V 



d 3 P 
(2tt) 3 



where S is the sum of the Boltzmann factors of the colour- 
singlet multiquark clusters allowed by the symmetries 

2N f 6N f -N c b 

S=J2 Cmb CX P H2me + N c b(e =p ft))0\ . (15) 

6—0 m— 5 fa 

In doing so we have replaced the ^ Q by the integra- 
tion over phase space volume d 3 xd 3 p/(2ir) 3 . It is in- 
teresting to note that due to the colour-singlet restric- 
tion the colour factor, N c = 3, and the baryon num- 
ber, 6, are always associated with the Boltzmann factor 
and thus with the partonic chemical potential indicating 
the excess number of quarks or antiquarks in addition 
to m number of quarks and antiquarks. This happens 
just because the diagonal matrix, R a a, in colour space 
or the Polyakov Loop re-adjusts the (anti) quarks of dif- 
ferent momenta to be in a same momentum state by al- 
lowing them to exchange the momentum in the thermal 
bath. As seen, in general, a given colour-singlet mode 
in (|15p has energy E m i, = (2m + N c b)e having parton 
content ](m + N c b)q, mq] for a hadron whereas that for 
anti-hadron is [(m+N c b)q, mq]. The pure mesonic modes 
(b = 0) can be identified with parton content (mq, mq) 
having total energy E m o = 2me. Similarly, the pure 
baryonic modes (m = 0) have the energy Eob = N c be 
with parton content [N c bq] for baryon whereas that for 
antibaryons is [N c bq]. In the above G m b is the weight 
factor, due to flavour and spin symmetry, appearing with 
each colour-singlet mesonic/baryonic/antibaryonic mode 
because of colour-singlet restriction. Their values for var- 
ious modes are listed in Ref. [l6l ]. For low lying mesons 
(to = 1 and 6 = 0) C\$ = 16 whereas for low lying 
baryons and antibaryons (to = and 6=1) Cox = 20, 
respectively. These are exact for the SU (2) flavor and the 
SU(2) spin symmetry [2(j of the quark model in which 
nucleons and deltas are degenerate. On the other hand 
to > 1 and 6 = correspond to excited mesonic modes 
whereas to > 1 and 6 > 1 are penta-quark and excited 
baryonic/ antibaryonic modes, which arc the Hagcdron 
states [2l| . 

Equation (|14p clearly exhibits a nontrivial result: 
SU(3) colour-singlet restriction, along with flavour 
and spin symmetry, on the quark-antiquark ensem- 
bles reorders the partition function in terms of 
Boltzmann factors of the colour-singlet multi-quark 
[mesonic/baryonic/antibaryonic] modes at any temper- 
ature. These are, however, not the bound states but 
can be regarded as aprecursor to the confinement due 
to recombination 0, 0, [HJ ■ Under a suitable confining 
mechanism (e.g., Polyakov Loop model [H, 13), one can 
hope that these multi-quark structures could evolve into 
colour-singlet hadrons in the low temperature limit. 
Probability and particle spectra: The probability of find- 
ing a single (anti)quark in the system in the energy in- 
terval e and e + de follows from (TT4"|) as [23| 



In [1 + S] 



(14) 



V(e)de 



Vd 3 p 



S] 



(16) 
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The logarithm is expanded with S < 1 that yields a so- 
lution e > (T, where ( ~ 1.7 for N j = 2. This provides 
an energy cut-off > 2£T for mesonic modes and > 3£T 
for baryonic modes. One can write the above as 



oo oo 



V(e)de=J2 J2 C *mb^P H 2 ™ e + 36(e T ft))0\ 

b—0 m— <5ob 
oo oo 

= E E KM* > 

6=0 m=(5ob 



(2^)3 
(17) 
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in which all higher order terms are accumulated in C* 
where C* nb = C m b for only low lying hadronic modes. In 
the above 7 ? ^ i) (e)de is the probability of a single parton, 
with energy e > £T, in the interval e and e + de in a given 
to6— th mode which reads as 



= ^nb exp [-(2me + 36(e T //))/?] 



(2tt) 3 ' 



(18) 



Now, the distribution of a parton within a given mb— th 
mode in the fluid in terms of the momentum of the parton 
{p^Un = e, is 4-velocity of the fluid) follows directly 
from (HU) as 



d 3 x d 3 p (2tt) 3 



exp [-(2me + 36(e =p m))/ 3 ] 



(19) 



One can easily invert the parton momentum distribution 



in the above into the momentum [P^u^ = E m b = 
36)e] distribution of a mb— th hadronic mode as 
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* e -(B m6 T3bM)/3 



(20) 



where n = (2m + 36) and the above is the Cooper-Frye 
distribution [24| at the freeze-out in the rest frame of 
the fluid. We assume that this distribution describes 
the corresponding hadrons after freeze-out. The above 
distribution has also an important consequence: the en- 
tropy would remain conserve since the number of quarks 
in quark's phase space is equal to that of hadrons in 
hadron's phase space. 

Now we can obtain the differential proton to pion ra- 
tio at central rapidity for a given transverse momentum 
P± » 3(T as 



dNi 



01 



dN 
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dN p 
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01 e 3 M /T = ± e HB/T 
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(21) 



where fis = 3/x, is the baryonic chemical potential and 
the above is independent of P±. PHENIX data show Q 
that the ratio has a plateau around unity in the P±_ range, 
2 GcV/c < P± < 4 GcV/c. Our estimation shows that 
it is in good agreement with RHIC data [H, [25|. Re- 
call that in the recombination model by Fries et al. @ 
it was found to be (5/3) exp(/is/T). Similarly, an- 
tiproton to pion is (5/4) exp(— hb/T) and that of an- 
tiproton to proton is exp(— 2/is/T), which are again in 



good agreement with the RHIC data [26|. Thus, the 
colour-singlet projection of the thermal parton ensemble 
based on symmetry consideration of the underlying the- 
ory also provides a strong basis for the hadronization by 
recombination of quark-gluon plasma at intermediate P± 
(2 GeV/c < Pj_ < 4 GeV/c) even though the ingredients 
are different from the recombination models 0, Q with a 
sudden approximation. 

Now the total number of a given hadron emitted by a 
fluid can be obtained from (TTS)) as 



N mb = 



VC* 



-nC±3bn/T 



1 + n( + 



n 2 C 2 



(22) 



which simply depends on the number of partonic degrees 
of freedom that forms the colour-singlet hadronic modes 
with e/T > (. The ratio of the proton to pion becomes 
~ 0.6 whereas it is ~ 0.3 for antiproton to pion with 
Hb/T ~ 0.33 in RHIC. This is consistent with the RHIC 
data [f|. 

Anisotropy and scaling: We now consider non-central col- 
lisions where the fluid has larger velocity on the x-axis 
(semi-minor) than on the y-axis (semi-major) leading to 
elliptic flow (2?| ■ The flow coefficient is defined as 



V2 



(cos 2(f>) 



(23) 



This requires that the distribution in (TT9")) should have a 
nontrivial dependence on azimuthal angle <f>. It is intro- 
duced through e = p^U/x = p±Uo(<f>) — p±u(<fi), where 
is the 4-velocity of the fluid and u(<f>) can be parameter- 
ized [13] in the following form 



u(<f>) = u + 2a cos 2(f> , 



(24) 



where u is 0— averaged of the maximum fluid 4— velocity 
in the <j) direction and a specifies the magnitude of the 
elliptic flow, which is about 4% in a non-central Au-Au 
collisions at RHIC. Considering u^u^ = Uq(</>) — w 2 (</>) = 
1 and expanding it to a first order in a, one can obtain 



uq(4>) = uq + 2av cos 2(f) 



(25) 



where v = u/uq. Now, the elliptic flow coefficient for the 
mb— th mode in (f2"3"|) is found to scale perfectly with the 
partonic p}_ as 



(«2) 



mb 



a , 
T 



v) (2m + 36) p ± 



(26) 



We note that this scaling is strictly valid in the re- 
combination region (2 GeV/c < P± < 4 GeV/c), where 
the mass dependence of the hadrons is irrelevant. How- 
ever, for P± < 1.5GeV/c the scaling deviates Q due to 
the mass ordering of the hadrons. To demonstrate this 
fact one needs a dynamical mass generating term of the 
quarks, which is, however, beyond scope of this calcula- 
tion. 

Summary: We propose that the recombination phenom- 
ena occurs readily when one constructs a colour-singlet 
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partition function from a thermal ensemble of quarks and 
antiquarks having internal symmetries, viz., colour, spin 
and flavour. This colour-singlet projection is shown to 
be related to the Polyakov Loop in the Polyakov gauge, 
which plays the important role for recombination. We 
also show that such recombination of thermal partons 
naturally describes some of the puzzling results in RHIC 
experiments, which, in turn, strongly supports that the 
recombination of thermal (anti)quarks could be one of 
the dominant mechanisms of hadronization in a dense 



medium in the intermediate transverse momentum of 
hadrons (2 GeV/c < P ± < 4 GeV/c). Nevertheless, to 
test the various quark recombination models for hadro- 
nisation from dynamical QCD still remain an open ques- 
tion. We, however, tried to address this question based 
on the symmetry consideration of the theory, which may 
be useful for the eventual solution of this important prob- 
lem. 
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